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Abstract. CaaCc^Og is a system exhibiting a series of fascinating properties, 
including magnetization plateaus and remarkably slow dynamics at low-T. These 
properties are intimately related to the geometrical frustration, which results from a 
particular combination of features: (i) the chains are arranged on a triangular lattice; 
(ii) there is a large uniaxial anisotropy; (iii) the intrachain and interchain couplings 
are ferromagnetic and antiferromagnetic, respectively. 

The uniaxial anisotropy is thus an issue of crucial importance for the analysis of the 
physical properties of Ca3Co206- However, it turns out that no precise investigation 
of this magnetic anisotropy has been performed so far. On the basis of susceptibility 
data directly recorded on single crystals, the present study reports on quantitative 
information about the anisotropy of CaaC^Og. 
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1. Introduction 

Ca3Co2C>6 is a geometrically frustrated spin-chain compound which has attracted 
considerable attention in recent years, owing to a series of puzzling properties 
PElElllElElEllHlEKinilllKlS]. The structure of this compound consists of chains 
made up of CoC>6 trigonal prisms alternating with C0O6 octahedra, which run along the 
c-axis of the hexagonal cell p] . These chains are separated by the Ca ions and they form 
an hexagonal lattice on the ab plane (see inset of Fig. [1]). The issue of the valence and 
spin states of the Co ions on these two sites has been the subject of intense controversies. 
Owing to a series of recent results p31 El E3 EEl Ej, it is now widely accepted that 
(i) the Co ions are trivalent for both sites; (ii) owing to the difference in the crystalline 
electric field (CEF), the Co 3+ (3<i 6 ) ions are in the high-spin state for the prismatic 
sites (S = 2) whereas they are in the low-spin state (S = 0) for the octahedral sites. 
In other respects, it was early found that the intrachain coupling J is ferromagnetic 
while the interchain coupling J' is antiferromagnetic [2]. Furthermore, CasQ^Oe was 
found to exhibit strongly anisotropic magnetic properties, the spins having a preferential 
orientation along the c axis [U [5] . Most probably, this last feature can be ascribed to 
the single-ion (uniaxial) anisotropy of the S = 2 spins at the prismatic sites. 

Ca 3 Co 2 6 possesses a combination of features (i.e. triangular lattice with spins 
oriented along the chain direction and with an antiferromagnetic interchain coupling) 
which yields a situation of geometrical frustration [IB] . It is clear that this frustration 
plays a great role in the peculiar magnetic properties of CasCo206. For instance, the 
long-range ordering which takes place below T/v — 26 K [2] was found to be closely 
related to the "Partially Disordered Antiferromagnetic" state (PDA) [3J proposed by 



Figure 1. Hysteresis loops (H || c) recorded with a sweep rate of 0.1 T / min. at 2 K 
(circles) and 10 K (solid line). The arrows indicate the direction of the field variation. 
The inset shows a projection of the structure along the hexagonal c-axis (the dark and 
light polyhedra represent Co06 trigonal prisms and CoC-6 octahedra, respectively; the 
grey circles represent the Ca ions; the solid lines emphasize the triangular arrangement 
of the chains in the ab plane). 
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Mekata [T5] about the ABX 3 compounds, another family of geometrically frustrated 
spin-chains [18] . In the PDA state, two over three chains are antiferromagnetically 
coupled (with long-range order along each of them), while the third one is left 
incoherent (i.e. without long-range order along its direction and no correlation with its 
neighbors). Among the peculiar magnetic properties of CasCc^Oe, one can mention (i) 
the appearance of a spin freezing at T S f (= ^/=o) ~ 8 K [3], O [7j , leading to a pronounced 
frequency dependence of the susceptibility (ATj/Tj)/Alog/ ~ 0.17) [12]; (ii) for 
T < Tsf, the existence of magnetization steps with a roughly constant field spacing 
[H [5j [20] and the appearance of a saturation in the spin-relaxation time, these features 
having led us to suggest the possibility of a phenomenon of Quantum Tunneling of the 
Magnetization (QTM) [20l [T2]. 

It can be emphasized that the geometrical frustration in Ca3Co206 is strongly 
related to the existence of a pronounced uniaxial anisotropy. In spite of this, we observe 
that no precise quantitative analysis of this magnetic anisotropy has been performed so 
far. It is the goal of the present study to extract quantitative information about the 
anisotropy of Ca 3 Co 2 6 . 

2. Experimental details 

Single crystals of CasCc^Og were grown according to the following method: a mixture 
of Ca3Co40g and K2CO3, in a weight ratio of 1/7, was heated at 950 °C for 50 hours 
in an alumina crucible in air; then, the cooling was performed in two steps, first 
down to 930 °C at 10 °C/h and then down to room temperature at 100 °C/h. This 
procedure leads to crystals having a needle-like shape (with the c axis along the longest 
dimension) which is convenient for orientation purposes. These crystals can be quite long 
(typically 4 mm) but they are thin (less than 0.5 mm). Accordingly, we performed the 
measurements on an assembly of crystals. Owing to their shape, we emphasize that the 
alignment between them, as well as with respect to the field direction (i.e. c axis either 
parallel or perpendicular to the field direction) can be precisely achieved. Moreover, the 
measurements for both orientations were carried out with the same mounting system 
— just by rotating the set of crystals inside the sample rod — in order to keep constant 
the residual background signal. 

Curves of magnetization as a function of T were recorded in 0.1 T, for both 
orientations, by using a Superconducting Quantum Interference Device magnetometer 
(MPMS, Quantum Design). We checked that, with 0.1 T, we are in the linear regime 
of the M(H) curves for both orientations and over the whole investigated temperature 
range, 50 < T < 300 K (i.e. T >> Tn). In what follows, one considers the susceptibility 
curves obtained by dividing the magnetization by the measuring field, leading to X\\(T) 
(i.e. with H || c) and x±(T) (i.e. with H _L c). Two observations support the reliability 
of these data: first, we found that the XllCO is i* 1 perfect agreement with that previously 
obtained using bigger crystals (of different shape); second, the powder-like data obtained 
by calculating Xp = 2x_i_)/3 is found to be superimposed on the susceptibility curve 
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registered on a ceramic (see Fig. [2]). 

3. Results and Discussions 

Figure [2] shows the X\\(T) an d X±(T) data — as well as the derived Xp(T) curve — along 
with data recorded on a ceramic sample and on bigger crystals with H || c. These 
data well exhibit the considerable magnetic anisotropy of this compound. To be more 
quantitative, however, one needs a proper modelization of the susceptibilities. It turns 
out that it is far from being obvious in the case of CasQ^Oe which gathers several 
particularities. 

As previously discussed, we can assume that the anisotropy of CasCk^Og is basically 
a single-ion feature, originating from a combination of the CEF and spin-orbit effects. 
We can also simply consider a Heisenberg form for the intrachain coupling, while the 
interchain coupling (much weaker) can be neglected in a first approximation. In other 
respects, we note that a spin chain compound would require in principle a specific 
treatment of the magnetic coupling in order to account for the enhanced fluctuation 
effects in ID systems. 

To the best of our knowledge, there is no model in the literature which takes into 
account all the features of CasCk^Og. In what follows, we will consider two modelizations 
which can be applicable — to some extent — to the case of CasQ^Og. We emphasize 
that each of these models only accounts for a part of the features of CasCk^Og. It 
also deserves to be noted that they correspond to different approaches that can be 
complementary in the analysis of magnetism. 
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Figure 2. Susceptibility data recorded on single crystals of Ca3Co206, with the 
magnetic field applied either along the c-axis (xpara referred to as x\\ m the text) or 
perpendicular to the c-axis {xperp referred to as x± m the text). Also shown is the 
corresponding "powder" data derived from the general expression Xp = (X||+ 2xj_)/3. 
The small filled triangles represent the data directly recorded on a ceramic sample. 
The small filled stars represent the data previously recorded on another set of crystals 
(different morphology with a larger size) with H \\ c. 
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3.1. First model 

We found in the literature a model which made a precise treatment of the on-site 
magnetism in a situation very close to ours. This work carried out by Parkin and 
Friend [21] dealt with iron intercalates of dichalcogenides, with Fe 2+ (i.e., 3d 6 like Co 3+ ) 
in a trigonally distorted anionic environment. It turns out that the cubic part of the 
CEF splits the 5 D ground state, and generates a T 5 triplet at lower energy, which 
can be represented by a fictitious orbital angular momentum L. On this basis, Parkin 
and Friend [21 J performed a rigorous analysis by considering the following single-site 
Hamiltonian: 

Hi = -k B 5Ll t - k B X Lj ■ Sj + ^(-Li + 2S 4 ) B . (1) 

The first term corresponds to the effect of the non-cubic part of the CEF, the second 
term is the spin-orbit (SO) coupling, and the last term is the basic form of the Zeeman 
energy. For our compound with such an Hamiltonian, 5 is positive while A is negative. 
There is no magnetic coupling in this model but we note that a spin-spin interaction 
can be added — if necessary — by using a mean-field (MF) approximation (even though 
we know such a MF treatment is poorly suited to ID systems). 

Using the eigenf unctions and associated energies, Parkin and Friend calculated the 
susceptibilities using the Van Vleck formula. Extending their approach, we took into 
account all the levels of the lower triplet (whereas only the two lowest are considered 
in Ref. [21J), but it actually makes almost no difference. The expressions found for the 
susceptibilities can be written as: 

= Nfi 2 B S(S + l) 1 
X|1 ~ 3k B T 

_ Nn 2 B S{S + 1) 1 
X± ~ W B T 

where 



h 
z 



(2) 



fi 



(3) 



/ll = 25 exp(-2s) + 9 exp(-s) + 1 + exp(s) + 

2s 2 5s 2 

9exp(2s) + 16exp( V p) + 4exp( h p) 

p p 



(4) 



fx = (-8/s) exp(-2s) - (4/s) exp(-s) + (4/s) exp(s) + (5) 

2s 2 5s 2 
(8/s) exp(2s) + 16exp( hp) + 4exp( h p) , 



2[exp(— 2s) + exp(— s) + 1 + exp(s) + exp(2s) + 

2s 2 5s 2 6s 2 

exp( \-p) + exp(— — hp) + l/2exp( h p)] 

p P V 



(6) 
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with s = X/T and p = —S/T. All others parameters have their usual meaning. 

It can be emphasized that the exact value of 5 has only a tiny influence on the 
final results (we used 5 = —100 K hereafter). Roughly speaking, one can consider 
that the main role of the non cubic part of the CEF is to split the T 5 triplet, which 
allows the SO to generate anisotropy. In other respects, it must be realized that A 
cannot really be regarded as a free parameter in these expressions, since the spin orbit 
parameter is supposed to be an intrinsic characteristic of each ion, not expected to vary 
significantly with the environment. In the case of Co 3+ , A must be close to -145 cm -1 
|22j. Accordingly, it deserves to be noted that there is almost no adjustable parameter 
in this model, apart from a possible correction related to J. 

Figure [3] shows the raw data, along with the curves calculated for two values of 
A. With the expected value A = —145 cm -1 , the calculated curves do not lie far from 
the experimental ones, they have the right shape in the high-T range, but there is a 
sizeable shift from the data. Moreover, the model cannot account for the existence of a 
minimum in l/x±{T), as it is visible in the data around 120 K. We also observed that 
adding a correction related to the intrachain coupling can yield a very nice fitting for 
X\\ (T) but not for x±{T)- Strikingly, it can be noticed that remarkable fittings of both 
curves can be obtained with A = —165 cm -1 . As previously discussed, however, such 
a significant departure from the free ion value of A is suspect. In conclusion, one must 
state that this model cannot well describe the magnetic anisotropy of Ca3Co20g, even 
though it well illustrates the leading role of the spin orbit coupling. 
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Figure 3. Experimental X\\(T) and X-l{T) curves, together with the calculations 
using the model of Parkin and Friend with A = —145 cm^ 1 (solid lines) and A = —165 
cm -1 (dashed lines). 
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3.2. Second model 

Let us now consider a very different approach, approximating the susceptibility in the 
form of series expansion. This general method discards the specificities of the system 
under investigation (in contrast to the previous model), but it was found to be powerful 
in many cases, even though it is by definition limited to the high-T part of the data 
(i.e. applicable only for kT larger than all the energy terms relevant to the system). 
Yosida |23j has derived such a series expansion which appears to be suitable to the 
case of Ca 3 Co20 6 , in the sense that it accounts for both a magnetic coupling J and an 
anisotropy parameter D. The starting point is an Hamiltonian of the form: 

Hi = -ksDSl, - 2k B J S.S,- + gfiaSiB , (7) 

where the index j refers to the first neighboring spins. In such a picture, the Zeeman 
term involves a Lande factor (which can be an adjustable parameter), while the uniaxial 
anisotropy appears as an effective term applied to the spin. Even though it can carry 
some ambiguities, this form of Hamiltonian is the most common to deal with uniaxial 
anisotropy, and it is very useful in practice. We underline that the S parameter of the 
previous model strongly differs from D in the above equation (this latter parameter 
incorporating both the CEF and SO effects). 

The formula derived by Yosida can be written as follows: 

X|| = y^[l + 2|], (8) 

x± = y^[i-f], (9) 

In these expressions, 

- C is the usual Curie constant 

c Ng*fj? B S(S + 1) 
3k B 

- 9 is the characteristic temperature associated to J (i.e. the Curie-Weiss 
temperature) 

2zJS(S+l) 

3 1 ' 

(here z is the number of nearest neighbors, i.e. 2 in the case of intrachain coupling) 

- Q is a characteristic temperature associated to D 



Q = D 



2S(S+1) 1 



(12) 



15 10_ 

Note that, in our case, 6 > and Q > 0. 

In Eqs (8-9), they are three free parameters which should be adjusted to best fit 
to both the an d XJ-CO data. Actually, one can also use combinations of X||C0 
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Figure 4. The main panel shows the l/xp(T) curve derived from the data (see 
text), which exhibits a linear behavior at high temperature (solid line). The inset 
shows the temperature dependence of the susceptibility ratio (xii/xx)) together with 
the expectations of the Yosida's model for several values of D: from bottom to top, 
D = 140 K (dotted line), D — 150 K (solid line) and D = 160 K (dashed line). 



and x±{T) allowing to isolate the influence of some of these parameters. For instance, 
it can be noted that the parameter D disappears in the powder formula, yielding the 
usual form 

X|| + 2x± C 

Xv = ^— = Y^e- (13) 

As shown in Fig. HI such a behavior is found to be well obeyed by the data. Fitting 
to the data in the range 200 — 300 K leads to g ~ 2.08 and J ~ 4.5 K. 
Furthermore, the ratio (yn/xj_) is expected to only depend on D 

Yll 1 + 2^ 

— = —^r (14) 

XL \-% 

Inset of Fig. H] shows that there is a restricted range of D values — around 150 K — 
which leads to a reasonable fitting in the high-T range. 

The X|| (T) and x±(T) curves, corresponding to g — 2.08, D = 150 K and J = 4.5 K, 
are reported in Fig. [5j The agreement with the data is far from being perfect, but it is 
reasonable if one considers that the same set of parameters must fit to two curves (not to 
mention that perfect fittings are obtained if one authorizes different sets of parameters 
for each curve). Probably, a part of the deviation from the data can be attributed to the 
fact that the temperatures of the fitting range are not much larger that D. Moreover, 
this model provides a rough description of the magnetic coupling since the ID character 
is only taken into account through the value of z (= 2). 
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Figure 5. Experimental X\\{T) and X±(T) curves, together with the calculations 
using the Yosida's model for the displayed values of parameters (according to the 
assumptions of this model, the calculations are limited to the high-T range). 

4. Conclusion 

The anisotropic susceptibility of the geometrically frustrated spin chain Ca3Co206 has 
been measured on single crystals. Two models of the literature have been used to 
analyze these data. None of them was found to be perfectly satisfying. However, the 
first one clearly pointed to the role of the SO in this anisotropy [21] . while the second one 
provided us with a reliable estimate of the effective anisotropy parameter D ~ 150 K 
[23] . This value is of importance to go further into the analysis of the peculiar magnetic 
properties exhibited by this compound. 

The present study also pointed to the issues that should be addressed in priority 
to better describe X\\(T) and x±{T) in CaaCc^Oe: (i) owing to the large values of D, 
it is necessary to develop a reliable model for T < D; (ii) since the intrachain coupling 
seems to have an impact even at large T (when considering a MF approach) , it would be 
valuable to consider a modelization accounting for the ID character of this compound. 
These developments are presently under way. 
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